Type-2 fuzzy neural network-based cooperative control method for wastewater treatment process

ABSTRACT

A cooperative fuzzy-neural control method is designed in this present invention. Due to the difficulty for cooperatively controlling the concentrations of the dissolved oxygen and nitrate nitrogen in wastewater treatment process, a cooperative fuzzy-neural control method is investigated. In this proposed method, firstly, a interval type-2 fuzzy neural network is employed to construct the cooperative fuzzy-neural controller. Secondly, a parameter cooperative strategy is proposed to cooperatively optimize the global and local parameters of the cooperative fuzzy-neural controller to meet the control requirements. This proposed cooperative fuzzy-neural control method can cooperatively control the concentrations of the dissolved oxygen and nitrate nitrogen in wastewater treatment process. The results illustrate that the proposed cooperative fuzzy-neural control method can achieve the high control accuracy and guarantee the normal operations of wastewater treatment process under the different operation conditions.

CROSS-REFERENCE TO RELATED PATENT APPLICATIONS

This application claims priority to Chinese Patent Application No. 201910899611.7, filed on Sep. 23, 2019, entitled “A type-2 fuzzy neural network-based cooperative control method for wastewater treatment process,” which is hereby incorporated by reference in its entirety.

TECHNICAL FIELD

In this present invention, a type-2 fuzzy neural network-based cooperative control (T2FNN-CC) method is designed to achieve the cooperative control of the dissolved oxygen (DO) and nitrate nitrogen (NO₃—N) concentrations in wastewater treatment process (WWTP). The concentration of DO and NO₃—N is the key control variables in WWTP, and play an important role for improving the effluent quality and energy consumption. It is feasible to implement the proposed T2FNN-CC in WWTP, and control the concentrations of the DO and NO₃—N cooperatively. This method can ensure the normal operation of WWTP under the different operation conditions. The technology of this present invention belongs to the field of water research and intelligent control.

BACKGROUND

The rapid growth of the global population and the rapid development of manufacturing industry accelerate the consumption of water resources, and produce a large amount of wastewater. In recent years, many countries have introduced a variety of environmental legislation to regulate the emission processes of wastewater and improve the emission standards to protect the environment. Wastewater treatment plants need to meet more requirements to reach the discharge standards. Meanwhile, the variable operation conditions is a great challenge to the stable operation of WWTP. Therefore, the research results of the present invention has the important significance for ensuring the stable operation of WWTP.

As the key control variables, the DO and NO₃—N concentrations can directly control the biochemical reaction processes in WWTP. In the wastewater treatment plants, the activated sludge process, which includes the ammonification reaction, nitrification reaction and denitrification reaction, is the most popular wastewater treatment method and the pollutants in wastewater can be removed by degrading the organic matter. The ammonification reaction can convert the organic nitrogen compounds into ammonia nitrogen, the nitration reaction can synthesize the ammonia nitrogen into NO₃—N, and the denitrification reaction can decompose the NO₃—N into nitrogen. Although these above reactions can achieve the effective removal of organic matter, it is difficult to control these reaction processes. For example, the nitration reaction is a kind of aerobic reaction and the sufficient DO concentration will promote the degradation of ammonia nitrogen in the nitration reaction. On the contrary, the denitrification reaction requires an anoxic environment and the NO₃—N concentration can be removed by the denitrification reaction under the hypoxic conditions. It is difficult to control both the nitration and denitrification reactions due to these conflicting conditions. The high DO concentration will affect the nitration and denitrification reactions and the NO₃—N concentration is a primary factor for the denitrification reaction. Therefore, it is urgent to control the DO and NO₃—N concentrations cooperatively. The type-2 fuzzy neural network can use the type-2 fuzzy rule to express the change of the process variables in WWTP and deal with the strong nonlinear and high uncertainty of WWTP. The proposed T2FNN-CC has strong robustness and high control precision, and use a parameter cooperative strategy to optimize the global and local parameters of type-2 fuzzy neural network to improve the control response speed and realize cooperative control of the DO and NO₃—N concentrations. The proposed T2FNN-CC can ensure the stable operation of WWTP under the different operation conditions, and has the practical application value.

In this present invention, a T2FNN-CC method is designed for controlling the DO and NO₃—N concentrations of WWTP, where the type-2 fuzzy neural network controller is used to control the controlled variables, and the parameter cooperative strategy is employed to update the global and local parameters of the type-2 fuzzy neural network controller.

SUMMARY

In this present invention, a T2FNN-CC method is designed. This proposed T2FNN-CC method can solve the control values of the DO and NO₃—N concentrations, and utilize the aeration value and internal backflow value to achieve the cooperative control of the DO and NO₃—N concentrations. The proposed T2FNN-CC method can ensure that the discharge standards are met, and improve the stability of WWTP under the different operation conditions.

The present invention adopts the following technical scheme and implementation steps:

A type-2 fuzzy neural network-based cooperative control method for controlling the dissolved oxygen (DO) and nitrate nitrogen (NO₃—N) concentrations of wastewater treatment process (WWTP), wherein aeration value and internal backflow value are used as control variables, the DO and NO₃—N concentrations are used as controlled variables, the method comprising the following steps:

(1) design a type-2 fuzzy neural network (T2FNN) for controlling the DO and NO₃—N concentrations, the T2FNN contains five-layers: an input layer, a membership layer, a rule layer, a consequent layer and an output layer, wherein:

the input layer contains 4 input neurons and an input vector is: X(t)=[x ₁(t),x ₂(t),x ₃(t),x ₄(t)]^(T)  (1) where X(t) is the input vector of T2FNN at time t, x₁(t) is an error between a set-point and a measured value of DO concentration at time t, x₂(t) is an error variation between the set-point and the measured value of DO concentration at time t, x₃(t) is an error between a set-point and a measured value of NO₃—N concentration at time t, x₄(t) is an error variation between the set-point and the measured value of NO₃—N concentration at time t, T represents a revolution of the matrix and vector;

the membership layer contains P membership neurons and a neuron represents an interval type-2 membership function:

$\begin{matrix} {{{\underline{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} \leq {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} > {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \end{matrix} \right.} & (2) \end{matrix}$ $\begin{matrix} {{{\overset{\_}{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \leq {{\underline{c}}_{ij}(t)}} \\ {1,} & {{{\underline{c}}_{ij}(t)} < {x_{i}(t)} < {{\overset{\_}{c}}_{ij}(t)}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \geq {{\overset{\_}{c}}_{ij}(t)}} \end{matrix} \right.} & (3) \end{matrix}$ where P is the total number of membership neurons, P=4M, M is the total number of rule neurons, 1<M≤15, m _(ij)(t) is a lower membership value of ith input to jth rule neuron at time t, m _(ij)(t) is an upper membership value of the ith input to the jth rule neuron at time t, 0<m _(j)(t)<m _(ij)(t)≤1, e is a natural constant and e=2.7183, c _(ij)(t) is a lower center of the ith input with the jth rule neuron at time t, c _(ij)(t) is an upper center of the ith input with the jth rule neuron at time t, 1<c _(ij)(t)<c _(ij)(t)<5, σ_(ij)(t) is standard deviation of the ith input with the jth rule neuron at time t, i is the number of inputs in the T2FNN and i=1, 2, 3, 4, j is the number of rule neurons and j=1, 2, . . . , M;

the rule layer contains M rule neurons, an output of each neuron is:

$\begin{matrix} {{F_{j}(t)} = \left\lbrack {{{\underline{f}}_{j}(t)},{{\overset{\_}{f}}_{j}(t)}} \right\rbrack} & (4) \end{matrix}$ $\begin{matrix} {{{{\underline{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\underline{m}}_{ij}(t)}}},{{{\overset{\_}{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\overset{\_}{m}}_{ij}(t)}}}} & (5) \end{matrix}$ where F_(j)(t) is firing strength of jth rule neuron at time t, f _(j)(t) is lower firing strength of the jth rule neuron at time t, f _(j)(t) is upper firing strength of the jth rule neuron at time t, 0<f _(j)(t)<f _(j)(t)≤1;

the consequent layer contains 4 consequent neurons, an output of each neuron is:

$\begin{matrix} {{{{{\underline{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\underline{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\underline{f}}_{j}(t)}}},{{{\overset{\_}{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\overset{\_}{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\overset{\_}{f}}_{j}(t)}}}}{{h_{j}^{k}(t)} = {{\sum\limits_{i = 1}^{4}{{w_{ij}^{k}(t)}{x_{i}(t)}}} + {b_{j}^{k}(t)}}}} & (6) \end{matrix}$ where y _(k)(t) is a lower output of consequent neuron with respect to kth output neuron at time t, y _(k)(t) is an upper output of consequent neuron with respect to the kth output neuron at time t, hk j(t) is jth consequent factor with respect to the kth output at time t, wk ij(t) is a consequent weight of ith input with respect to jth rule neuron and the kth output at time t, bk j (t) is the deviation of the jth rule neuron with respect to the kth output at time t, k=1, 2;

the output layer contains 2 neurons, an output of each neuron is: Δu _(k)(t)−q ^(k)(t) y _(k)(t)+(1−q ^(k)(t)) y _(k)(t)  (7) where Δu_(k)(t) is an output of kth output neuron at time t, q^(k)(t) is a scale factor of the kth output neuron, 0<q^(k)(t)<1;

(2) train T2FNN, which includes:

1) divide parameters of T2FNN into global and local parameters, and define a global parameter vector and a local parameter vector as: Φ_(g)(t)=[ c _(ij)(t), c _(ij)(t),σ_(ij)(t)] Φ_(l)(t)=[w _(ij) ^(k)(t),b _(j) ^(k)(t),q ^(k)(t)]  (8) where Φ_(g)(t) is the global parameter vector at time t, Φ_(l)(t) is the local parameter vector at time t, and an objective function is:

$\begin{matrix} {{{l(t)} = {{{\alpha(t)}{l^{1}(t)}} + {\left( {1 - {\alpha(t)}} \right){l^{2}(t)}}}}{{\alpha(t)} = \frac{❘{l^{1}(t)}❘}{{❘{l^{1}(t)}❘} + {❘{l^{2}(t)}❘}}}{{l^{1}(t)} = {{y_{d}^{1}(t)} - {y^{1}(t)}}}{{l^{2}(t)} = {{y_{d}^{2}(t)} - {y^{2}(t)}}}} & (9) \end{matrix}$ where l(t) is a composite error at time t, α(t) is an error coefficient at time t, l¹(t) is the error between the set-point and measured value of DO concentration at time t, l²(t) the error between the set-point and measured value of NO₃—N concentration at time t, y1 d(t) is the set-point of DO concentration at time t, y2 d(t) is the set-point of NO₃—N concentration at time t, y¹(t) is a real value of DO concentration at time t, y²(t) is a real value of NO₃—N concentration at time t;

2) utilize parameter cooperative strategy to optimize the global and local parameters of T2FNN cooperatively, an updating rule of parameters is:

$\begin{matrix} {{{\Phi\left( {t + 1} \right)} = {{\Phi(t)} + {\left( {{H(t)} + {{\theta(t)}I}} \right)^{- 1}{G_{l}(t)}}}}{{H(t)} = {{J^{T}(t)}{J(t)}}}{{G_{l}(t)} = {{J^{T}(t)}{l(t)}}}{{\theta(t)} = {{❘\frac{l(t)}{{l(t)} + {l\left( {t - 1} \right)}}❘}{\theta\left( {t - 1} \right)}}}{{J\left( \text{⁠}t \right)} = {\left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{\Phi_{g}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\Phi_{l}(t)}}} \right\rbrack = \left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{{\underline{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{{\overset{\_}{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\sigma_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{w_{ij}^{k}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{b_{ij}^{k}(t)}}\text{⁠},\text{⁠}\frac{\partial{l(t)}}{\partial{q^{k}(t)}}} \right\rbrack}}} & (10) \end{matrix}$ where Φ(t+1) is a parameter vector at time t+1, Φ(t) is a parameter vector at time t, I is a unit matrix, H(t) is a quasi Hessian matrix at time t, G_(l)(t) is an error gradient vector at time t, θ(t) is an adaptive learning rate at time t and θ(t)∈(0, 1], J(t) is a Jacobian vector at time t, ∂l(t)/∂Φ_(g)(t) is a partial derivative of the composite error with respect to the global parameter vector at time t, ∂l(t)/∂Φ_(l)(t) is a partial derivative of the composite error with respect to the local parameter vector at time t, ∂l(t)/∂c _(ij)(t) is partial derivatives of the composite error with respect to lower uncertain center at time t, ∂l(t)/∂c _(ij)(t) is partial derivatives of the composite error with respect to the upper uncertain center at time t, ∂l(t)/∂σ_(ij)(t) is partial derivatives of the composite error with respect to the standard deviation at time t, ∂l(t)/∂wk ij(t) is partial derivative of the composite error with respect to the consequent weight at time t, ∂l(t)/∂bk j(t) is partial derivative of the composite error with respect to the deviation at time t, ∂l(t)/∂q^(k)(t) is partial derivative of the composite error with respect to the scale factor at time t;

(3) design the type-2 fuzzy neural network-based cooperative control method for controlling the DO and NO₃—N concentrations of WWTP, which includes:

1) calculate outputs of T2FNN according to Eq. (7);

2) compare a value of composite error with a pre-set threshold, if l(t)>β, go to step 3), and if l(t)≤β, go to step 4), where β is the pre-set threshold and β=0.01;

3) calculate updated values of parameters according to Eq. (10);

4) calculate control outputs at current time: u _(k)(t)=u _(k)(t−1)+Δu _(k)(t)  (11) where u_(k)(t) is kth control output at time t, u₁(t) is the control output with respect to an aeration value at time t, u₂(t) is the control output with respect to an internal backflow value at time t, u_(k)(t−1) is kth control output at time t−1,

5) enter u_(k)(t) into WWTP to control the DO and NO₃—N concentrations, return to step 1);

(4) use u₁(t) and u₂(t) to control the concentration of DO and NO₃—N in WWTP, u₁(t) is the control input of the aeration value at time t, u₂(t) is the control input of the internal backflow value at time t, the control results are the concentrations of DO and NO₃—N in WWTP.

The Novelties of this Present Disclosure Contain:

(1) WWTP is a complex industrial process with the strong nonlinearity, high uncertainties and multiple disturbances, and the exact mathematic model of WWTP is difficult to design according to the general methods. Since the artificial neural network owns the strong learning ability and approximation capability, the T2FNN-CC method can be employed to achieve the precise and stable control for WWTP under the different operation conditions.

(2) In this present invention, a T2FNN-CC method is studied to realize the cooperative control of the DO and NO₃—N concentrations. This proposed T2FNN-CC method can utilize the robustness of T2FNN and the parameter cooperative strategy to meet the control requirements for WWTP under the different operations conditions, and address the difficulty of the stability control in WWTP.

Attention: for the convenience of the description, the proposed T2FNN-CC method is used to control the DO and NO₃—N concentrations of WWTP cooperatively. All of the control methods, which utilize the principle of this present invention, also belong to the scope of this present invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows the control scheme of this present invention.

FIG. 2 shows the structure of T2FNN in T2FNN-CC method

FIG. 3 shows the control results of DO concentration in T2FNN-CC method.

FIG. 4 shows the control errors of DO concentration in T2FNN-CC method.

FIG. 5 shows the control results of NO₃—N concentration in T2FNN-CC method.

FIG. 6 shows the control errors of NO₃—N concentration in T2FNN-CC method.

DETAILED DESCRIPTION

In this present invention, a T2FNN-CC method is designed. This proposed T2FNN-CC method can solve the control values of the DO and NO₃—N concentrations, and utilize the aeration value and internal backflow value to achieve the cooperative control of the DO and NO₃—N concentrations. The proposed T2FNN-CC method can ensure that the discharge standards are met, and improve the stability of WWTP under the different operation conditions.

The present invention adopts the following technical scheme and implementation steps:

A type-2 fuzzy neural network-based cooperative control method for controlling the dissolved oxygen (DO) and nitrate nitrogen (NO₃—N) concentrations of wastewater treatment process (WWTP), wherein aeration value and internal backflow value are used as control variables, the DO and NO₃—N concentrations are used as controlled variables, the control scheme is shown in FIG. 1 , the method comprising the following steps:

(1) design a type-2 fuzzy neural network (T2FNN) for controlling the DO and NO₃—N concentrations, the T2FNN contains five-layers: an input layer, a membership layer, a rule layer, a consequent layer and an output layer, the structure of T2FNN is shown in FIG. 2 and the details of the layers:

wherein:

the input layer contains 4 input neurons and an input vector is: X(t)=[x ₁(t),x ₂(t),x ₃(t),x ₄(t)]^(T)  (1) where X(t) is the input vector of T2FNN at time t, x₁(t) is an error between a set-point and a measured value of DO concentration at time t, x₂(t) is an error variation between the set-point and the measured value of DO concentration at time t, x₃(t) is an error between a set-point and a measured value of NO₃—N concentration at time t, x₄(t) is an error variation between the set-point and the measured value of NO₃—N concentration at time t, T represents a revolution of the matrix and vector;

the membership layer contains P membership neurons and a neuron represents an interval type-2 membership function:

$\begin{matrix} {{{\underline{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} \leq {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} > {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \end{matrix} \right.} & (2) \end{matrix}$ $\begin{matrix} {{{\overset{\_}{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \leq {{\underline{c}}_{ij}(t)}} \\ {1,} & {{{\underline{c}}_{ij}(t)} < {x_{i}(t)} < {{\overset{\_}{c}}_{ij}(t)}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \geq {{\overset{\_}{c}}_{ij}(t)}} \end{matrix} \right.} & (3) \end{matrix}$ where P is the total number of membership neurons, P=4M, M is the total number of rule neurons, 1<M≤15, m _(ij)(t) is a lower membership value of ith input to jth rule neuron at time t, m _(ij)(t) is an upper membership value of the ith input to the jth rule neuron at time t, 0<m _(ij)(t)<m _(ij)(t)≤1, e is a natural constant and e=2.7183, c _(ij)(t) is a lower center of the ith input with the jth rule neuron at time t, c _(ij)(t) is an upper center of the ith input with the jth rule neuron at time t, 1<c _(ij)(t)<c _(ij)(t)<5, σ_(ij)(t) is standard deviation of the ith input with the jth rule neuron at time t, i is the number of inputs in the T2FNN and i=1, 2, 3, 4, j is the number of rule neurons and j=1, 2, . . . , M;

the rule layer contains M rule neurons, an output of each neuron is:

$\begin{matrix} {{F_{j}(t)} = \left\lbrack {{{\underline{f}}_{j}(t)},{{\overset{\_}{f}}_{j}(t)}} \right\rbrack} & (4) \end{matrix}$ $\begin{matrix} {{{{\underline{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\underline{m}}_{ij}(t)}}},{{{\overset{\_}{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\overset{\_}{m}}_{ij}(t)}}}} & (5) \end{matrix}$ where F_(j)(t) is firing strength of jth rule neuron at time t, f _(j)(t) is lower firing strength of the jth rule neuron at time t, f _(j)(t) is upper firing strength of the jth rule neuron at time t, 0<f _(j)(t)<f _(j)(t)≤1;

the consequent layer contains 4 consequent neurons, an output of each neuron is:

$\begin{matrix} {{{{{\underline{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\underline{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\underline{f}}_{j}(t)}}},{{{\overset{\_}{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\overset{\_}{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\overset{\_}{f}}_{j}(t)}}}}{{h_{j}^{k}(t)} = {{\sum\limits_{i = 1}^{4}{{w_{ij}^{k}(t)}{x_{i}(t)}}} + {b_{j}^{k}(t)}}}} & (6) \end{matrix}$ where y _(k)(t) is a lower output of consequent neuron with respect to kth output neuron at time t, y _(k)(t) is an upper output of consequent neuron with respect to the kth output neuron at time t, hk j(t) is jth consequent factor with respect to the kth output at time t, wk ij(t) is a consequent weight of ith input with respect to jth rule neuron and the kth output at time t, bk j (t) is the deviation of the jth rule neuron with respect to the kth output at time t, k=1, 2;

the output layer contains 2 neurons, an output of each neuron is: Δu _(k)(t)=q ^(k)(t) y _(k)(t)+(1−q ^(k)(t)) y _(k)(t)  (7) where Δu_(k)(t) is an output of kth output neuron at time t, q^(k)(t) is a scale factor of the kth output neuron, 0<q^(k)(t)<1;

(2) train T2FNN, which includes:

1) divide parameters of T2FNN into global and local parameters, and define a global parameter vector and a local parameter vector as: Φ_(g)(t)=[ c _(ij)(t), c _(ij)(t),σ_(ij)(t)] Φ_(l)(t)=[w _(ij) ^(k)(t),b _(j) ^(k)(t),q ^(k)(t)]  (8) where Φ_(g)(t) is the global parameter vector at time t, Φ_(l)(t) is the local parameter vector at time t, and an objective function is:

$\begin{matrix} {{{l(t)} = {{{\alpha(t)}{l^{1}(t)}} + {\left( {1 - {\alpha(t)}} \right){l^{2}(t)}}}}{{\alpha(t)} = \frac{❘{l^{1}(t)}❘}{{❘{l^{1}(t)}❘} + {❘{l^{2}(t)}❘}}}{{l^{1}(t)} = {{y_{d}^{1}(t)} - {y^{1}(t)}}}{{l^{2}(t)} = {{y_{d}^{2}(t)} - {y^{2}(t)}}}} & (9) \end{matrix}$ where l(t) is a composite error at time t, α(t) is an error coefficient at time t, l¹(t) is the error between the set-point and measured value of DO concentration at time t, l²(t) the error between the set-point and measured value of NO₃—N concentration at time t, y1 d(t) is the set-point of DO concentration at time t, y2 d(t) is the set-point of NO₃—N concentration at time t, y¹(t) is a real value of DO concentration at time t, y²(t) is a real value of NO₃—N concentration at time t;

2) utilize parameter cooperative strategy to optimize the global and local parameters of T2FNN cooperatively, an updating rule of parameters is:

$\begin{matrix} {{{\Phi\left( {t + 1} \right)} = {{\Phi(t)} + {\left( {{H(t)} + {{\theta(t)}I}} \right)^{- 1}{G_{l}(t)}}}}{{H(t)} = {{J^{T}(t)}{J(t)}}}{{G_{l}(t)} = {{J^{T}(t)}{l(t)}}}{}{{\theta(t)} = {{❘\frac{l(t)}{{l(t)} + {l\left( {t - 1} \right)}}❘}{\theta({t - 1})}}}{{J\left( \text{⁠}t \right)} = {\left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{\Phi_{g}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\Phi_{l}(t)}}} \right\rbrack = \left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{{\underline{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{{\overset{\_}{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\sigma_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{w_{ij}^{k}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{b_{ij}^{k}(t)}}\text{⁠},\text{⁠}\frac{\partial{l(t)}}{\partial{q^{k}(t)}}} \right\rbrack}}} & (10) \end{matrix}$ where Φ(t+1) is a parameter vector at time t+1, Φ(t) is a parameter vector at time t, I is a unit matrix, H(t) is a quasi Hessian matrix at time t, G_(l)(t) is an error gradient vector at time t, θ(t) is an adaptive learning rate at time t and θ(t)∈(0, 1], J(t) is a Jacobian vector at time t, ∂l(t)/∂Φ_(g)(t) is a partial derivative of the composite error with respect to the global parameter vector at time t, ∂l(t)/∂Φ_(l)(t) is a partial derivative of the composite error with respect to the local parameter vector at time t, ∂l(t)/∂c _(ij)(t) is partial derivatives of the composite error with respect to lower uncertain center at time t, ∂l(t)/∂c _(ij)(t) is partial derivatives of the composite error with respect to the upper uncertain center at time t, ∂l(t)/∂σ_(ij)(t) is partial derivatives of the composite error with respect to the standard deviation at time t, ∂l(t)/∂wk ij(t) is partial derivative of the composite error with respect to the consequent weight at time t, ∂l(t)/∂bk j(t) is partial derivative of the composite error with respect to the deviation at time t, ∂l(t)/∂q^(k)(t) is partial derivative of the composite error with respect to the scale factor at time t;

(3) design the type-2 fuzzy neural network-based cooperative control method for controlling the DO and NO₃—N concentrations of WWTP, which includes:

1) calculate outputs of T2FNN according to Eq. (7);

2) compare a value of composite error with a pre-set threshold, if l(t)>β, go to step 3), and if l(t)≤β, go to step 4), where β is the pre-set threshold and β=0.01;

3) calculate updated values of parameters according to Eq. (10);

4) calculate control outputs at current time: u _(k)(t)=u _(k)(t−1)+Δu _(k)(t)  (11) where u_(k)(t) is kth control output at time t, u₁(t) is the control output with respect to an aeration value at time t, u₂(t) is the control output with respect to an internal backflow value at time t, u_(k)(t−1) is kth control output at time t−1,

5) enter uk(t) into WWTP to control the DO and NO3—N concentrations, return to step 1);

(4) use u1(t) and u2(t) to control the concentration of DO and NO3—N in WWTP, u1(t) is the control input of the aeration value at time t, u2(t) is the control input of the internal backflow value at time t, the control results are the concentrations of DO and NO3—N in WWTP. FIG. 3 gives the control results of DO concentration, X axis shows the time, and the unit is day, Y axis is the DO concentration, and the unit is mg/L, the black solid line is the set-points of DO concentration, and the black dotted line is the real value of DO concentration. FIG. 4 gives the control errors of DO concentration, X axis shows the time, and the unit is day, Y axis is the errors of DO concentration, and the unit is mg/L. FIG. 5 gives the control results of NO₃—N concentration, X axis shows the time, and the unit is day, Y axis is the NO₃—N concentration, and the unit is mg/L, the black solid line is the set-points of N₃—N concentration, and the black dotted line is the real value of NO₃—N concentration. FIG. 6 gives the control errors of NO₃—N concentration, X axis shows the time, and the unit is day, Y axis is the errors of NO₃—N concentration, and the unit is mg/L. 

What is claimed is:
 1. A type-2 fuzzy neural network-based cooperative control method for controlling dissolved oxygen (DO) and nitrate nitrogen (NO₃—N) concentrations of wastewater treatment process (WWTP), wherein aeration value and internal backflow value are used as control variables, the DO and NO₃—N concentrations are used as controlled variables, the method comprising the following steps: (1) design a type-2 fuzzy neural network (T2FNN) for controlling the DO and NO₃—N concentrations, the T2FNN contains five-layers: an input layer, a membership layer, a rule layer, a consequent layer and an output layer, wherein: the input layer contains 4 input neurons and an input vector is: X(t)=[x₁(t), x₂(t), x₃(t), x₄(t)]^(T)  (1) where X(t) is the input vector of the T2FNN at time t, x₁(t) is an error between a set-point and a measured value of DO concentration at time t, x₂(t) is an error variation between the set-point and the measured value of DO concentration at time t, x₃(t) is an error between a set-point and a measured value of NO₃—N concentration at time t, x₄(t) is an error variation between the set-point and the measured value of NO₃—N concentration at time t, T represents a revolution of the weight matrix of the input layer and the input vector; the membership layer contains P membership neurons and a neuron represents an interval type-2 membership function: $\begin{matrix} {{{\underline{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} \leq {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},{{x_{i}(t)} > {\left( {{{\underline{c}}_{ij}(t)} + {{\overset{\_}{c}}_{ij}(t)}} \right)/2}}} \end{matrix} \right.} & (2) \end{matrix}$ $\begin{matrix} {{{\overset{\_}{m}}_{ij}\left( {x_{i}(t)} \right)} = \left\{ \begin{matrix} {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\underline{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \leq {{\underline{c}}_{ij}(t)}} \\ {1,} & {{{\underline{c}}_{ij}(t)} < {x_{i}(t)} < {{\overset{\_}{c}}_{ij}(t)}} \\ {e^{{- \frac{1}{2}}{(\frac{{x_{i}(t)} - {{\overset{\_}{c}}_{ij}(t)}}{\sigma_{ij}(t)})}^{2}},} & {{x_{i}(t)} \geq {{\overset{\_}{c}}_{ij}(t)}} \end{matrix} \right.} & (3) \end{matrix}$ where P=4M, M is a number of rule neurons, 1<M≤15, m_(y)(t) is a lower membership value of ith input to jth rule neuron at time t, m _(ij)(t) is an upper membership value of the ith input to the jth rule neuron at time t, 0<m _(ij)(t)<m _(ij)(t)≤1, e is a natural constant and e=2.7183, c_(ij)(t) is a lower center of the ith input with the jth rule neuron at time t, c _(ij)(t) is an upper center of the ith input with the jth rule neuron at time t, 1<c _(ij)(t)<c _(ij)(t)<5, σ_(ij)(t) is standard deviation of the ith input with the jth rule neuron at time t, i is the number of inputs in the T2FNN and i=1, 2, 3, 4, j is the number of rule neurons and j=1, 2, . . . , M; the rule layer contains M rule neurons, an output of each neuron is: $\begin{matrix} {{F_{j}(t)} = \left\lbrack {{{\underline{f}}_{j}(t)},{{\overset{\_}{f}}_{j}(t)}} \right\rbrack} & (4) \end{matrix}$ $\begin{matrix} {{{{\underline{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\underline{m}}_{ij}(t)}}},{{{\overset{\_}{f}}_{j}(t)} = {\prod\limits_{i = 1}^{4}{{\overset{\_}{m}}_{ij}(t)}}}} & (5) \end{matrix}$ where F_(j)(t) is firing strength of jth rule neuron at time t,f _(j)(t) is lower firing strength of the jth rule neuron at time t, f _(j)(t) is upper firing strength of the jth rule neuron at time t, 0<f _(j)(t)<f _(j)(t)≤1; the consequent layer contains 4 consequent neurons, an output of each neuron is: $\begin{matrix} {{{{{\underline{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\underline{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\underline{f}}_{j}(t)}}},{{{\overset{\_}{y}}_{k}(t)} = \frac{\sum\limits_{j = 1}^{M}{{{\overset{\_}{f}}_{j}(t)}{h_{j}^{k}(t)}}}{\sum\limits_{j = 1}^{M}{{\overset{\_}{f}}_{j}(t)}}}}{{h_{j}^{k}(t)} = {{\sum\limits_{i = 1}^{4}{{w_{ij}^{k}(t)}{x_{i}(t)}}} + {b_{j}^{k}(t)}}}} & (6) \end{matrix}$ where y _(k)(t) is a lower output of consequent neuron with respect to kth output neuron at time t, y _(k)(t) is an upper output of consequent neuron with respect to the kth output neuron at time t, hk j(t) is jth consequent factor with respect to the kth output at time t, wk ij(t) is a consequent weight of ith input with respect to jth rule neuron and the kth output at time t, bk j (t) is the deviation of the jth rule neuron with respect to the kth output at time t, k=1, 2; the output layer contains 2 neurons, an output of each neuron is: Δu _(k)(t)=q ^(k)(t) y _(k)(t)+(1−q ^(k)(t)) y _(k)(t)  (7) where Δu_(k)(t) is an output of kth output neuron at time t, q^(k)(t) is a scale factor of the kth output neuron, 0<q^(k)(t)<1; (2) train the T2FNN, which includes: 1) divide parameters of the T2FNN into global and local parameters, and define a global parameter vector and a local parameter vector as: Φ_(g)(t)=[ c _(ij)(t), c _(ij)(t),σ_(ij)(t)] Φ_(l)(t)=[w _(ij) ^(k)(t),b _(j) ^(k)(t),q ^(k)(t)]  (8) where Φ_(g)(t) is the global parameter vector at time t, Φ_(l)(t) is the local parameter vector at time t, and an objective function is: $\begin{matrix} {{{l(t)} = {{{\alpha(t)}{l^{1}(t)}} + {\left( {1 - {\alpha(t)}} \right){l^{2}(t)}}}}{{\alpha(t)} = \frac{❘{l^{1}(t)}❘}{{❘{l^{1}(t)}❘} + {❘{l^{2}(t)}❘}}}{{l^{1}(t)} = {{y_{d}^{1}(t)} - {y^{1}(t)}}}{{l^{2}(t)} = {{y_{d}^{2}(t)} - {y^{2}(t)}}}} & (9) \end{matrix}$ where l(t) is a composite error at time t, α(t) is an error coefficient at time t, l¹(t) is the error between the set-point and measured value of DO concentration at time t, l²(t) the error between the set-point and measured value of NO₃—N concentration at time t, yl d(t) is the set-point of DO concentration at time t, y2 d(t) is the set-point of NO₃—N concentration at time t, y¹(t) is a real value of DO concentration at time t, y²(t) is a real value of NO₃—N concentration at time t; 2) utilize parameter cooperative strategy to optimize the global and local parameters of the T2FNN cooperatively, an updating rule of parameters is: $\begin{matrix} {{{\Phi\left( {t + 1} \right)} = {{\Phi(t)} + {\left( {{H(t)} + {{\theta(t)}I}} \right)^{- 1}{G_{l}(t)}}}}{{H(t)} = {{J^{T}(t)}{J(t)}}}{{G_{l}(t)} = {{J^{T}(t)}{l(t)}}}{{\theta(t)} = {{❘\frac{l(t)}{{l(t)} + {l\left( {t - 1} \right)}}❘}{\theta\left( {t - 1} \right)}}}{{J\left( \text{⁠}t \right)} = {\left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{\Phi_{g}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\Phi_{l}(t)}}} \right\rbrack = \left\lbrack \text{⁠}{\frac{\partial{l(t)}}{\partial{{\underline{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{{\overset{\_}{c}}_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{\sigma_{ij}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{w_{ij}^{k}(t)}},\text{⁠}\frac{\partial{l(t)}}{\partial{b_{ij}^{k}(t)}}\text{⁠},\text{⁠}\frac{\partial{l(t)}}{\partial{q^{k}(t)}}} \right\rbrack}}} & (10) \end{matrix}$ where Φ(t+1) is a parameter vector at time t+1, Φ(t) is a parameter vector at time t, I is a unit matrix, H(t) is a quasi Hessian matrix at time t, G_(t)(t) is an error gradient vector at time t, θ(t) is an adaptive learning rate at time t and θ(t)∈(0, 1], J(t) is a Jacobian vector at time t, ∂l(t)/∂Φ_(g)(t) is a partial derivative of the composite error with respect to the global parameter vector at time t, ∂l(t)/∂Φ_(l)(t) is a partial derivative of the composite error with respect to the local parameter vector at time t, ∂l(t)/∂c ^(ij)(t) is partial derivatives of the composite error with respect to lower uncertain center at time t, ∂l(t)/∂c _(ij)(t) is partial derivatives of the composite error with respect to the upper uncertain center at time t, ∂l(t)/∂σ_(ij)(t) is partial derivatives of the composite error with respect to the standard deviation at time t, ∂l(t)/∂wk ij(t) is partial derivative of the composite error with respect to the consequent weight at time t, ∂l(t)/∂bk j(t) is partial derivative of the composite error with respect to the deviation at time t, ∂l(t)/∂q^(k)(t) is partial derivative of the composite error with respect to the scale factor at time t; (3) design the type-2 fuzzy neural network-based cooperative control method for controlling the DO and NO₃—N concentrations of the WWTP, which includes: 1) calculate outputs of the T2FNN according to Eq. (7); 2) compare a value of composite error with a pre-set threshold, if l(t)>β, go to step 3), and if l(t)≤β, go to step 4), where β is the pre-set threshold and β=0.01; 3) calculate updated values of parameters according to Eq. (10); 4) calculate control outputs at current time: u _(k)(t)=u _(k)(t−1)+Δu _(k)(t)  (11) where u_(k)(t) is kth control output at time t, u₁(t) is the control output with respect to an aeration value at time t, u₂(t) is the control output with respect to an internal backflow value at time t, u_(k)(t−1) is kth control output at time t−1, 5) enter u_(k)(t) into the WWTP to control the DO and NO₃—N concentrations, return to step 1); (4) use u₁(t) and u₂(t) to control the concentration of DO and NO₃—N in the WWTP, u₁(t) is the control input of the aeration value at time t, u₂(t) is the control input of the internal backflow value at time t, the control results are the concentrations of DO and NO₃—N in the WWTP. 